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Abstract — We analyze the feasibility of distributed spatial 
multiplexing with limited feedback in a slow-fading interference 
network with MN non-cooperating single-antenna sources and 
M non-cooperating single-antenna destinations. In particular, we 
assume that the sources are divided into M mutually exclusive 
groups of N sources each, every group is dedicated to transmit a 
common message to a unique destination, all transmissions occur 
concurrently and in the same frequency band and a dedicated 
1-bit broadcast feedback channel from each destination to its 
corresponding group of sources exists. We provide a feedback- 
based iterative distributed (multi-user) beamforming algorithm, 
which "learns" the channels between each group of sources 
and its assigned destination. This algorithm is a straightforward 
generalization, to the multi-user case, of the feedback-based 
iterative distributed beamforming algorithm proposed recently by 
Mudumbai, Hespanha, Madhow and Barriac in [1] for networks 
with a single group of sources and a single destination. Putting 
the algorithm into a Markov chain context, we provide a simple 
convergence proof. We then show that, for M finite and N — > oo, 
spatial multiplexing based on the beamforming weights produced 
by the algorithm achieves full spatial multiplexing gain of M 
and full per-stream array gain of N, provided the time spent 
"learning" the channels scales linearly in N. The network is 
furthermore shown to "crystallize" in the sense that, in the large- 
N limit, the M individual fading links not only decouple (as 
reflected by full spatial multiplexing gain) but also converge 
to non-fading links. Finally, we quantify the impact of the 
performance of the iterative distributed beamforming algorithm 
on the crystallization rate, and we show that the multi-user 
nature of the network leads to a significant reduction in the 
crystallization rate, when compared to the M = 1 case. 

I. INTRODUCTION 

We consider a special class of interference networks, where 
MN non-cooperating sources are divided into M mutually 
exclusive groups Qui = 1, 2, . . . , M, such that the N sources 
in the i t h group, denoted as Sf, j = 1, 2, . . . , N, are dedicated 
to transmit, through slow-fading channels, a common message 
to their assigned single-antenna destination Dj. All transmis- 
sions occur concurrently and in the same frequency band and 
the destinations T>i do not cooperate. This network models the 
second hop of the coherent multi-user relaying protocol in [2] 
under the assumption that the first hop transmission is error- 
free. The results in [2] imply that, for the interference network 
considered in this paper, for M fixed and N — > oo, full spatial 
multiplexing gain of M and a per-stream (distributed) array 
gain of N can be obtained, provided that each source knows 
the channel to its assigned destination perfectly. In this paper, 
we analyze the case where the perfect channel state informati- 
on assumption is relaxed to having a 1-bit broadcast feedback 
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channel from each destination T>i to its sources Qi. These 
broadcast feedback channels are non-interfering. We provide a 
feedback-based iterative distributed (multi-user) beamforming 
algorithm, which "learns" the channels between each group 
of sources and its assigned destination. This algorithm is a 
straightforward generalization, to the multi-user case, of the 
feedback-based iterative distributed beamforming algorithm 
proposed recently by Mudumbai, Hespanha, Madhow and 
Barriac in [1] for networks with a single group of sources 
and a single destination. Making the simplifying assumption, 
compared to [1], of the fading coefficients as well as all 
the signals being real-valued allows us to put the iterative 
algorithm into a Markov chain context, thereby setting the 
stage for a simple convergence proof. We then show that, 
for M finite and N — > oo, spatial multiplexing based on 
the beamforming weights produced by the iterative algorithm 
achieves full spatial multiplexing gain of M and full per- 
stream (distributed) array gain of N, provided the time spent 
"learning" the channels scales linearly in N. We furthermore 
demonstrate that the M effective links Qi — > T>i in the network 
not only decouple (reflected by full spatial multiplexing gain) 
but also converge to non-fading links as TV — > oo, i.e., in 
the terminology of [3], the network "crystallizes". Finally, we 
quantify the impact of the performance of the iterative algo- 
rithm on the crystallization rate, and we show that the multi- 
user nature of the network leads to a significant reduction in 
the crystallization rate, when compared to the M = 1 case. 

Notation: The superscripts T and _1 stand for transposition 
and inverse, respectively. \Q\ denotes the cardinality of the set 
Q, \x\ is the absolute value of the scalar x, and [a\ denotes 
the greatest integer that is smaller than or equal to the real 
number a. A/"(^i,cr 2 ) stands for the normal distribution with 
mean fi and variance a 2 , log(-) denotes logarithm to the base 
2. fx(-) stands for the probability density function (p.d.f.) of 
the random variable X and X <~ Y denotes equivalence in 
distribution. A = B denotes that the sets (of terminals) A and 
B are equal. P(u>) is the probability of event u, E[X] and 
VAR[X] are the expected value and the variance, respectively, 
of the random variable X and w.p. stands for with probability. 
Since the terminals in Qi are assumed to have a common 
message for their assigned destination T>i, we will be using 
the notation Qi — > Dj to denote the corresponding single- 
input single-output link between the group Qi and destination 
T>i. Vectors and matrices are set in lower-case and upper-case 
bold-face letters, respectively. 

II. System and Signal Model 

We assume that Xi[n] is the common message of the sources 
in Qi to be transmitted to T>i. In the remainder of the paper, we 
distinguish between a training phase during which the N scalar 



channels between S(,j = 1,2, ... ,N, and 2?, are "learned" 
for each i £ {1,2,..., M}, and a data transmission phase fol- 
lowing the training phase. During the training phase for Qi the 
feedback broadcast channel T). L — > Qi is used once every frame 
of Tf time slots. The l t h frame (I = 0, 1, ...), denoted by Tu 
consists of the time slots n = ITf, ITf + 1, . . . , (1 + 1)T/ - 1. 
In the Z t /j frame, each source Sj multiplies the sequence Xi [n] 
(which is a training sequence during the training phase) by a 
corresponding beamforming weight a\[l] £ {1,-1} before 
transmission; these beamforming weights are kept constant 
during the entire frame. We furthermore assume that all the 
channels in the network are flat-fading and remain constant 
throughout the entire time-interval of interest, i.e., during 
training and data transmission phases. Throughout the paper, 
we assume that M is finite. The corresponding input-output 
relations are now given by 
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in the n t h time slot, h\ r € 



denotes the symbol received at T>i 

R stands for the fading coefficient 
between S 3 and T>i and Wi[n] denotes the 7V(0, N a ) i.i.d. noise 
sequence at T>i. The fading coefficients h\ r ,for all i,r,j, are 
assumed i.i.d. Af(0, 1). The signals Xi[n) obey the average 
power constraint 
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so that the average power transmitted by each group Qi 
is limited by P. During the training phase, the frame-rate 
sequences a\ [l] are updated based on the 1-bit feedback 
received at the end of each frame. The goal of this process 
is to find the beamforming weights a 3 = sign(ft,^ i ). While it 
is in general difficult to put the individual groups into perfect 
beamforming configuration, so that 
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we will show that if in each group a sufficient number of 
sources is in beamforming configuration, full spatial multi- 
plexing gain can be achieved. In what follows, we shall say 
that in a given group Qi, sources satisfying h 3 i o? i > are 
aligned, whereas sources with h\ i p? i < are reverse-aligned. 
We denote the final beamforming weights a\ produced during 
the training phase as ol\. In the data transmission phase, 
the groups Gi,i — 1,2, ...,M, perform co-channel data 
transmission (i.e., spatial multiplexing) based on eel so mat 
the corresponding input-output relation is given by ([T]i with 
oA \l] = oA . 



III. Iterative Distributed Beamforming 

We shall next describe the algorithm carried out during the 
training phase. The overall training phase is assumed to consist 
of T tr — k a MN frames (each of which contains Tf time slots) 
divided into M blocks of k a N frames each. The role of the 
parameter k Q will become clear later. During each of these M 
blocks, precisely one of the groups Qi follows the three-step 
iterative distributed beamforming algorithm, described below, 
while all the other groups of sources Q r ,r ^ i, remain silent. 
At the end of the training phase of k MN frames, each 
of the groups Qi is in (close-to) beamforming configuration 
with respect to (w.r.t.) its assigned destination. The order in 
which the groups follow the three-step procedure below can 
be decided offline and communicated to all the nodes in the 
network. Without loss of generality (w.l.o.g.), we assume that 
the group Qi is being processed during the i t h block defined as 
the set of frames I = (i—l)k N, (i—l)k N+l, . . . , ik N—l. 
Since much of the analysis in the current and the next section 
deals with a single group Qi only, we shall consider, w.l.o.g., 
the group Q\, drop the index i = 1 and, wherever appropriate, 
use the convention Q = Q\ , = <Sf , T> = D\ , ot 1 = o? y , a? = 
a{ and h? =h\ x . The three steps carried out by the iterative 
distributed beamforming algorithm can now be summarized as 
follows. 

• Step 1. Initialization of the received signal level: This 
step pertains to the zeroth frame in the block. Each of the 
sources S 3 initializes its beamforming weight according 
to a J [0] = 1, initializes an auxiliary beamforming weight 
as a? [0] = ct' J [0] and starts transmitting the pilot symbol 
yj P/N . The corresponding received signal at destination 
T> is given by 
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The destination T> then estimates the received signal level 
by averaging y[n] over the entire frame, resulting in 
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Here, we assumed that the estimate of L rx is perfect, 
which requires that Tf be sufficiently large. Finally, V 
initializes L mdx = U x . 

Step 2. Iterative distributed beamforming: This is the 
iterative step that is performed for each frame, except 
for the zeroth one, i.e., for I = 1,2,..., k Q N — 1 (recall 
that the zeroth frame is used to carry out the previous 
initialization step). The details of this step are as follows. 
At the beginning of each frame, each of the sources 
S 3 ,j = 1,2,..., AT, chooses its beamforming weight 



a J [l], independently across j, according to: 
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Each of the sources S 3 then transmits the pilot symbol 
J P/N throughout the frame, using the beamforming 
weight a J [l], i.e., S J transmits the signal a 3 [l]yP~/N . At 
the end of the frame under consideration, T> estimates the 
corresponding received signal level, as in the initialization 
step, according tcQ 



(3) 



and, through the 1-bit broadcast feedback channel, in- 
forms all the sources in Q whether L IX > L mm or not. 
Based on the received feedback, the sources 6> J , j = 
1, 2, . . . , TV, update their auxiliary beamforming weights 




' [I] as follows: 



a 3 \l] = 



\& 3 \l-l], if U x < L mdx 

• ~ , j = 1,2,. . . , N. 

\a 3 [l], if L a > L max 



(4) 

Finally, if U x > L mdx , V updates L mix = U x . 
• Step 3. Silencing: At the end of the block, the sources 
in Q store the current values of their respective auxiliary 
beamforming weights as a 3 — a 3 [k N — 1] and go silent. 
In the proposed protocol, only one group of sources is 
active during a given block. It is therefore natural to ask 
whether the groups that have finished their training phase could 
start their data transmission phase while the remaining groups 
are "learning" their channels through the three-step procedure 
above. At first sight, one would be tempted to conclude 
that such a modified protocol would result in higher spectral 
efficiency. We note, however, that performing Steps 1 and 2 
above in the presence of interference created by the groups 
already transmitting data would require a longer Tf, say Tj, 
in order to achieve the same quality of the received signal 
level estimate as in the original protocol, i.e., in the absence of 
interference. The resulting tradeoff can be illustrated roughly 
by assuming maximum-likelihood estimation of the parameter 
U x . The corresponding variance of the estimation error is 
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in the original protocol and by a 
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Consequently, if we require the variance of the estimation error 
in the two protocols to be equal, then 
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Assuming that in each of the links already in (close-to) be- 
amforming configuration the corresponding group of sources 
transmits at a rate of R bits per time slot, the total number 
of bits transmitted in the modified protocol during the first 
k MNTj time slots is given by Rk Q NTjM(M - l)/2 [bit] 
whereas it is Rk M 2 N(Tj -Tf) [bit] in the original protocol. 
Using |6]), this implies that the modified protocol would have 
a higher spectral efficiency if 



< 1 - 
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Now, since the weights a\ have been optimized to be in 
beamforming configuration w.r.t. h\ p they are independent 
of h\ r for r ^ i. Furthermore, a\ = sign(/i| i ) and 
ol\ = — siga(hl 4 ) for 5| in the set of aligned and reverse- 
aligned sources, respectively. Since P(siga(h J ii ) = 1) = 
P(sign(/i| i ) = -1) = 1/2, it follows that P{5P i = 1) = 
P(oi = —1) = 1/2, for all j. Along with the assumption 



h\ r ~ Af(0, 1), for all i, r, j, we therefore get h\ r o? r ~ h\ r , 
for all 7 and r 

2 
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< i, which allows us to conclude that 

,&l , r < i, scales linearly in N. Consequently, it 

follows from (5j that a] is proportional to P, which implies 
that the condition Q would not be met for any reasonable 
choice of SNR = P/N Q . We can therefore conclude that, 
provided the variance of the signal level estimation error is the 
relevant performance measure, the original protocol would in 
practice always outperform the modified protocol in terms of 
spectral efficiency. 

IV. Convergence of The Iterative Distributed 
Beamforming Algorithm 

We shall next show that the iterative algorithm described 
in the previous section converges to the beamforming confi- 
guration for k a large. The proof is rather straightforward and 
consists of the following sequence of arguments: 

We start by recognizing that through (|2]l and Q, the co- 



efficient vector a[l] = [a 1 ^] 



a N [l}} T depends 



only on the coefficient vector in the previous frame, i.e., 
on at[l— 1]. The sequence of vectors a[l] therefore follows 
a Markov chain. 

« Each of the elements in the vector a[l] can take on the 
values 1 or —1, which implies that the total number of 
states of the Markov chain is given by 2^. Let us denote 
these states by Qi, Qz,- ■■ , <22 N - To each of the states, 
we can associate a value YljLi hP&P . 

• The Markov chain has one special state, namely, when 
oP — sign(/i J ), for all j. W.l.o.g. we let this state be 
Qi and note that it corresponds to the beamforming 
configuration, i.e., 



'Again, we assume the estimate to be perfect. 
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which implies that once the system is in state Qi the 
parameter L max is maximized over all states Qi,i — 
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I, 2,..., 2 N . From the update rule Q, we can therefore 
conclude that once in Qi the Markov chain will remain 
in Qi, which implies that Qi is an absorbing state (see 
[4, Def. 11.1]). 

• Each of the states Q r ,r ^ 1, corresponds to a coeffi- 
cient vector a with at least one reverse-aligned element. 
Let us denote the number of reverse-aligned elements 
corresponding to Q r by r a , where r a > 1. Then, the 
transition probability from Q r to the absorbing state Qi 
is (1//V) r °(l - l/N) N - r °. We can therefore conclude 
that the absorbing state can be reached from all the states 
Q r so that the Markov chain is an absorbing Markov 
chain (see [4, Def. 11.1]). It then follows from [4, Th. 

II. 3] that the system moves to the absorbing state Q\ 
w.p. 1 as the length of the training phase becomes large. 
Since Q\ corresponds to the beamforming configuration, 
we can conclude that the iterative algorithm eventually 
converges to the beamforming configuration. 

We emphasize that due to our simplifying assumptions, 
compared to [1], of the fading coefficients and the signals 
being real-valued, we were able to prove convergence of 
the iterative algorithm to the beamforming configuration in a 
simple fashion using basic results from Markov chain theory. 
The proof above, however, does not reveal anything about 
the rate of convergence of the algorithm. It seems difficult 
to obtain results on the actual convergence rate because the 
state-transition probabilities in the Markov chain, implicitly, 
depend on the actual realizations of the fading coefficients 
h? . However, interesting insights into the behavior of the 
convergence rate, as a function of N, can be obtained by 
considering convergence in expectation according to 
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In particular, we shall show that convergence in the sense of |9) 
can be obtained if the length of the training phase T tr scales 
linearly in N, i.e., k is independent of N. To be precise, we 
shall establish convergence to the beamforming configuration 
up to a certain level in the sense that we will allow for 
a (small) fraction of the sources to be reverse-aligned. The 
corresponding concept and the associated convergence proof 
are provided next. 

A. Convergence to e-level and convergence proof 

In the following, we shall be interested in convergence of 
the iterative algorithm in the sense of |9]) up to a certain level. 
Concretely, we shall allow that, on average, a (small) fraction 
e of the sources in a given group is reverse-aligned, which 
trivially implies that (1 — e) is the fraction of sources that are 
aligned, on average. Throughout the paper, we shall assume 
that e is independent of N ancQeiV € Z + . Denoting the sets 

2 Strictly speaking, for a given N, this requires that e be an integer multiple 
ofl/JV. 



of aligned and reverse-aligned sources after the update Q at 
the end of the t t h frame as A[t] and A[t], respectively, we 
have 

5>'o>[f]= e \ hj \- E i^i- ( 10 > 

53 eg sieA[t] s^eA[t] 

We say that convergence in the sense of |9]) up to an e-level 
has been achieved if 
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We can now formalize our convergence result as follows. 
Theorem 1: For any e G > and large TV, the number 
of iterations required in the second step of the distributed 



beamforming algorithm in Section III to achieve convergence 
to e -level is at most k N, where k a is a constant independent 
of N. 

Proof: We start by defining 

s new M = E - E WdP ^ - 1 \ (13) 

sjeg sieg 
so that S" lew [i] > 0, for all t. With this definition, we have 
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Let there be (l-e)TV and eN sources in ^4[t— 1] and!A[t-l], 
respectively, with < e Q < e < 1. Moreover, let q a [t] and 
q a [t] denote the number of sources in A[t — 1] and A[t — 1], 
respectively, that alter their beamforming weights from 1 to 
— 1 or vice-versa according to |2| at the beginning of the tth 
frame. Then, we get 
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We next show that the expected value inside the summation 
in ( fT~6] > can be lower-bounded by c a s, where c Q is a constant 
independent of N. To this end, we start by noting that q a [t] = 
0, q a [t] = s corresponds to the case where precisely s sources 



move from A[t — 1] to A[t] and none of the sources moves 
from A[t — 1] to A[t]. We denote the sets of these s sources 
as 7Z, and note that since the s sources in 7Z are chosen 
from the eN sources in A[t — 1], there are precisely ( e ^) 
possible choices for 1Z, with the corresponding sets denoted 
as 7£i,7\!.2, . . . ,TZ/eN\. Furthermore, each source in A[t — 1] 
alters its beamforming weight independently and with equal 
probability, and hence, each of the ( £ ^) choices is equally 
likely. From ( fT3) l it therefore follows that for a given set TZi, 
we have 
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where Step (a) is a result of the fact that each source in — 1] 
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of the sets TZi. We therefore get 
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where Step (a) follows from \A[t — 1]| = eN and the fact 
that the h? are identically distributed. Next, we need to show 
that the integral on the right hand side (RHS) of fj"8] > can be 
lower-bounded by a constant independent of N. To this end, 
we start by noting that 

= e f\M\\Xlt-l]( X ' 
which, using e > e G , implies 



/|W|0) 



(19) 



We can then argue that finding the minimum of J xf(x) dx 
over the class of functions /(x) that satisfy 

/(#) = 0, x < 0, 
f(x) > 0, x > 0, 
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guarantees that this minimum is a lower bound on 
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class of functions. Concretely, we want to determine 
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where the minimization is under the constraints ( |20"] >. This 
minimization problem can be solved as follows. We start by 
setting x = Q _1 (^^ 11 ), where Q{v) = ^=e~~^ du, 

so that using f\hi\{x) = 
0, x < 0, we have 
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Since the second term in the last line of ( |2"2"| i is positive due 
to the constraint f(x) < (l/e )f\f l j\(x), it follows that setting 
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0, otherwise 

yields the desired lower bound. Further, corresponding to f(x) 
in ( |23] l, we get 
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Substituting ( |24] i into ( fl"8j ) and the result thereof into ( fT6| ), we 
obtain 



[ E W[/]}* =0 |W [<S' new [t]]] 

eN 

>c o P(q a [t]=0)Y / sP{<lalt]=s) 

1 



iV 



s=0 
(l-e)JV 



E- 

.9 = 1 



eJV 



1-1 
N 



1 

N eN 



1 - 



tN-s 



= Co 1 ~ 



iV 



E' 

8=1 

iV eJV 



N 



I"* 
1 

at - i 



AT 

if 



JV 



1 



JV - 1 



eA 7 1 + 



1 



JV — 1 



siV-1 



(25) 
(26) 
(27) 

(28) 
(29) 
(30) 



(since V ( K ) kx^ 1 = K(l + x)^ 1 ) 
k=i W 



AT 



AT/ AT- 1 V N - 1 



:JV-1 



ft, e 1 



> c D e e" 



1 

iV 

-(l-£) 
"(l-^o) 



(l-e)JV 



(for large N) 
(since e > e Q ). 



(31) 

(32) 

(33) 
(34) 



Substituting ( |34| ) into ( fTB) finally yields 

k N-l 

> T — 



E 



{" J l']}i°o \ hJ 



E* 



eg 



1 = 1 



(l-e„) 



which upon setting 



(l-2e )e 1 - e ° 



-E[|^'| 



and noting that k Q does not depend on A~ establishes the 
desired result. □ 

B. Numerical Results 

Next, we present simulation results related to the con- 
vergence behavior of the iterative distributed beamforming 
algorithm described in Section [TTT] In particular, Fig[T] shows 
how E[^ 5jeg /i J a J [t]l (averaged over 50 realizations of the 
fading coefficients /V) evolves as a function of t for A~ = 
100, 500, 1000 and 5000, respectively. In all four cases, con- 
vergence occurs within approximately 10 A 7 iterations, thereby 
corroborating the fact that the convergence time is linear in 
AT. 



v. achievability of multiplexing gain and 
Crystallization 

The aim of this section is to show that performing data 
transmission, i.e., spatial multiplexing, based on the beam- 
forming weights 5i\ obtained from the iterative distributed 



beamforming algorithm described in Section III results in a 



spatial multiplexing gain of M and a per-stream (distributed) 
array gain proportional to N. Moreover, using the framework 
introduced in [3], we prove that the network "crystallizes", i.e., 
the individual fading links Qi — > Dj converge to non-fading 
links as A" — > oo. Finally, we quantify the impact of M and 
e D on the crystallization rate, i.e., the rate, as a function of 
N, at which the fading links Q. L — > T>i converge to non-fading 
links. 

Consider the group Qi and let Ai and Ai denote the set 
of aligned and reverse-aligned sources, respectively, corre- 
sponding to the beamforming weights a\, so that \A%\ = 
(1 — e a )N and \Ai\ = e N. Then, the corresponding input- 
output relations for the individual links Qi — * T>i are 



.</.» •'•.'»: E EK 

M j N v 

+ E (EHr^J^W +Wi[n], 1 = 1,2, 

r=l ^ j=l ' 

r^i 



,M. 



interference 

Note that here we assume that for any given realization of the 
fading coefficients h\ i7 precisely (1 — eo)AT of the sources are 
aligned. Convergence of the iterative distributed beamforming 
algorithm in expectation according to \\2) , however, guaran- 
tees only that this is the case on average. We furthermore 
assume Gaussian codebooks in what follows. Under the ass- 
umption of the destination T>. L having perfect knowledge of its 
effective channel coefficient (J^. Wa\ ~ I Mi I) an< ^ °f 



the coefficients Xw=i h 3 ir a J r ,for all r ^ i, corresponding to 
the effective interference channels, the outage probability for 
the link Qi — > T>i is given by 



/'( '/.o-:.! • SIMM • R 



P (SINRj < 2 



= P 




< 2 2R > - 1 



< 2 



211.. 



(36) 
(37) 

(38) 
(39) 



To upper-bound the outage probability, we use the (union) 
bounding techniques introduced in [3]. Skipping the details, 
we note that for any positive k\, k 2 and fc 3 such that k\ > fc 2 , 
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Fig. 1 . Convergence of the iterative distributed beamforming algorithm. 



the following upper bound holds: 



p , J^J^I-£i>mIJ < ( fcl - k 2 y 

<p\ [ E i^i -E i^i) <( fc i-M 2 

2 



N N 

+ >fc 3 



E p 



E%< 



> 



M - 1 



-4. ' Ai 



E^ c 



> 



M - 1 



(41) 



(42) 



(40) 



where we implicitly assumed that M > 1 and &3 > (NN ) /P. 
Next, we use the fact, noted previously at the end of Secti- 
on ITT, that h\ r a? r ~ /i^ r , for r ^ i. Consequently, we have 

l fe lJ-Ex I^mD 2 < (fci - fc 2 ; 



V I V h j aP? 4- ^ 

<p(j2\ h U< k ^)+ p (j2\ h U>^ 



A, 



+ (M-1)P 



> 



M - 1 



(43) 



Each of the three terms in ( |43| l can be upper-bounded using 
large deviations bounds [5], In particular, employing 



P (EN > fc ) - ( 2 
P ^ fc ) ; 



!e 



in (|43]l, we get 



pr(^) < 



efci 



(l-Co)W 



(l-e )iVV ^ 
+ 2(M - l) e "2«(M-i) . 



+ 2e 2 <^«) 2 



(44) 



The key to obtaining meaningful upper bounds from |44) 
lies in a judicious choice of the constants k\ , k 2 and fc 3 
ensuring that, for Pj = (l/2)log(l + (fci - k 2 ) 2 /k 3 ), 
P° ut (Ri) — > as A — > 00 while satisfying the conditions 

fci > fe 2 and fc 3 > (NN a )/P. Since e[£x Kill « A 



cx A, motivated by the above large 



and VAR 

deviations bounds, it is sensible to set 

(l-e )(N-VN) pit 



(45) 



Similarly, since Y1~J IM il deviates around a mean value 
proportional to A with a variance proportional to TV and 
I J2j h\ r I deviates around a mean value proportional to y/~N 
with a variance proportional to A, again motivated by the 
above large deviations bounds, we set 



V / 2(l + ln2)(e ^ + Vn) 



and k 3 — (M - 1)A 



1+6 



NN 
P 



(46) 



with the constant 5 > 0. Note that the condition k\ > k% 



implies that £q < 



0.2419. With the above 



l+eV^( 1+ln2 ) 

choices for the parameters hi, k 2 and ^3, in the limit TV — > 00, 
we get 

o2A-l _ - fc 2) 2 



so that Pj = (1/2) log (1 + cxiV 1 " 15 ). Substituting (|45j» and 
( |46] i into ( |44) , in the large- limit, we finally obtain 



P out [ -log(l + Cl A 



1-^ 



< 1 



+ e 



- eo iV-2^ + 2(M _ 1)e -M 



e -VW + e -e oN -2VlV + 2(M _ 1)e -^- 



(48) 



We can therefore conclude that P° M (Ri) — > as A — > 00 
for any rate P, < (l/2)log(l + c x A) (recall that J can be 
arbitrarily small). Since this holds true for all groups Q i7 we 
can choose Ri = (l/2)log(l + CjJV) with a < c\, for all i, 
and get P° ut (P,) — > 0, for all i, as A — > 00, which implies 
full spatial multiplexing gain of M, a per-stream array gain 
proportional to A, and convergence of each of the links Qi — > 
I?i to a non-fading link. In summary, in the language of [3], 
we can conclude that the network "crystallizes" as A — > 00. 
The third term on the RHS of ( |48] i nicely reflects the impact 
of interference on the crystallization rate. Specifically, for 6 < 
1/2, this term dominates the decay rate as a function of A. A 
smaller S corresponds, through P; = (1/2) log(l + c\ A 1- * 5 ), 
to higher data rates, but results in a reduced crystallization 
rate. In the single-user case, i.e., for M = 1, the third term 
equals zero reflecting the absence of interference. We can 
therefore conclude that the crystallization rate in the presence 
of interference is significantly smaller than in the single-user 
case M = 1. Finally, regarding the proportionality constant 
c%, it can be observed that the smaller e (i.e., the smaller the 
fraction of reverse-aligned sources) the larger c\, and hence 
the larger the individual rates P; = (l/2)log(l + qA) still 
guaranteeing crystallization. On the other hand, for smaller 
e a the second term on the RHS of (|48]l becomes larger, again 
reflecting that a higher data rate comes at the cost of increased 
outage probability. 
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